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Scattering experiments have led to many important 
discoveries in physics. In the field of dilute atomic gases, 
an understanding of inter-atomic scattering, as well as 
scattering between light and atoms, is fundamental to 
our ability to prepare, control, and manipulate systems 
for study. For example, a Feshbach resonance is a scat- 
tering resonance that makes it possible to control the 
strength and sign of inter-atomic interactions pQ. An 
understanding of this phenomenon has led to the forma- 
tion of molecular condensates [21 [3] and the investigation 
of the crossover between Fermi superfluidity and Bose- 
Einstein condensation [4j [5] . 

In this paper we consider elastic single-particle scatter- 
ing from a one-dimensional two-component degenerate 
Fermi gas when the projectile particle is identical to one 
of the confined species. Our theoretical treatment is a 
number-conserving approach based on the Hartree-Fock 
ground state of the trapped gas. We construct the excita- 
tions in a configuration-interaction approach using one- 
particle continuum excitations and bound two-particle 
one-hole excitations. The excitations are single-particlc- 
like in the asymptotic limit and by extracting the scat- 
tering phase shifts we determine the scattering properties 
of the system. 

We predict Fano-type scattering resonances that arise 
due to the interrelation of the one- and multi-particle 
branches of the excitation spectrum. In particular, when 
the energy of a two-particle one-hole bound state lies 
within the one-particle continuum, coupling between the 
two branches significantly modifies the scattering proper- 
ties of the system near the uncoupled bound state energy. 
In this system, the two branches of the excitation spec- 
trum are coupled by the inter-atomic pairing. 

We also present a BCS mean-field description of the 
scattering. In this approach, the scattering properties of 
the system are determined from the asymptotic behav- 
ior of the many-body quasiparticle excitations. We find 
that, in the uncoupled system, bound hole excitations 
lie within the particle excitation continuum. The mean- 
field pair potential couples the particle and hole branches 
of the excitation spectrum and we again observe Fano- 
type scattering resonances. However, these resonances 



are quantitatively different from those predicted by the 
configuration-interaction method. In this particular sys- 
tem, the BCS treatment leads to spurious results that 
are explicitly linked with the violation of particle num- 
ber conservation in the theory. Effects of this type are 
particularly evident in small systems. 

Our paper is organized as follows. In Sec. [T] we out- 
line the system for study. In Sec. |H]we investigate elas- 
tic single-particle scattering from a trapped Fermi gas 



III 



using the configuration-interaction method. In Sec 
we present the BCS mean-field treatment of that same 
problem. In Sec. |IV| we give a detailed discussion of the 
validity of the mean-field approach, and we conclude in 
Sec.lVl 



I. TWO-COMPONENT FERMI GAS 

We consider a trapped one-dimensional degenerate 
Fermi gas at zero temperature with two equally popu- 
lated spin components interacting via an attractive con- 
tact potential. The Hamiltonian is 



(i) 



where the field operator if> a (x) destroys a particle at po- 
sition x, in the spin state a =T,4, and g < 0. The 
single-particle Hamiltonian is 



(2) 



where M is the atomic mass and V ext (x) is a symmet- 
ric external trapping potential with a zero energy con- 
tinuum threshold. The exact form of the trapping po- 
tential is not crucial for our discussion. However, for 
our numerical calculations we use the Gaussian form 
Vext(x) = — Vb exp(— 2x 2 /w 2 ), with Vq > 0, which can 
be realized for atoms in a waveguide by applying a sin- 
gle Gaussian laser beam. The Gaussian width provides 
a convenient energy scale E w — ti 2 /2Mw 2 . 
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II. CONFIGURATION-INTERACTION 
METHOD 

To describe elastic single-particle scattering we deter- 
mine the excitation spectrum using the equation of mo- 
tion method [B]. The ground state |G) of the trapped 
gas has energy Eq, i.e., H\G) = Eq\G). We introduce 
the operators Qt that create excitations \v) = Ql\G) 
satisfying H\v) = E y \u), and the excitation spectrum is 
given by 

[ff,Qt]|G) = ^Qt|G), (3) 

where E v = E v — Eq. Projecting Eq. ^ onto the state 
A"t|G), where X is an arbitrary operator, gives 

(G\X[H,Ql]\G) = E v (G\XQl\G). (4) 

Later, X is chosen so that Eq. Q generates a matrix 
eigenvalue equation for the excited states \v). Equation 
Q is exact but is difficult to solve because we do not 
know the ground state |G) or the excitation creation op- 
erators Ql,. 

An approximate solution to Eq. Q can be found by 
neglecting pair correlations in the ground state, i.e., we 
calculate the matrix elements using the Hartree-Fock 
ground state |HF) in place of |G). We briefly summa- 
rize the Hartree-Fock method here. The Hartree-Fock 
Hamiltonian is 

Hhf = J2 [ ^(x) [H SP {x) + W{x)] Mx)dx, (5) 

where W(x) — g(HF\ipl l (x)'ip a (x) |HF) is spin indepen- 
dent because we have chosen equal spin populations. The 
expansion ip a (x) = ~Y^ n <f> n {x)a na diagonalizes Hamilto- 
nian ([5]) where the Hartree-Fock wavefunctions satisfy 

[H SP (x) + W(x)} = E n <t> n (x). (6) 

The operators cv na and a na obey fermionic commutation 
relations and the modes arc populated according to 

(BF\dl ia a n 'p\BF) = n„<W<W, (7) 

where, at zero temperature, n n = 1 for E n < E-p and 
n n = for E n > Ep. We choose the Fermi energy Ep < 
so that all of the particles are bound. The Hartree-Fock 
ground state is constructed by adding one particle of each 
spin to the lowest available energy level until the Fermi 
energy is reached, i.e., 

I HF >= I II ^Al) i°>> ( g ) 

where np denotes the highest occupied level. 

We determine the Hartree-Fock ground state of the 
system by numerically calculating the self-consistent 
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FIG. 1: (a) The Hartree-Fock ground state of a trapped Fermi 
gas. The curves correspond to (thick) the Hartree poten- 
tial W(x), (thin) the Gaussian external potential V ex t{x), and 
(dashed) the combined potential V ax t{x) + W(x). (b) A scat- 
tering one-particle excitation, and (c) a bound two-particle 
one-hole excitation, where the energy level occupation of the 
(•) spin up and (o) spin down particles is indicated schemati- 
cally. The arrows indicate the particle rearrangement between 
the configurations in panels (b) and (c). (d) The ground state 
energy level occupation [see Eq. |7|] is indicated by the hori- 
zontal component of the markers x . In all panels, the horizon- 
tal lines indicate (dashed) the Fermi energy E-p and (solid) the 
bound Hartree-Fock energies E n < [see Eq. H]. Parameters 
are g = -9.55™£„, V = 127.32S„, and E F = -19.10E W . 

Hartree potential W(x). We use an iterative method 
where the kinetic energy term in Eq. Q is evaluated ac- 
cording to a finite-difference formula. Figure [lja) shows 
the Hartree potential for a particular set of parame- 
ters where there are six particles of each spin state, i.e., 
{N a ) = J vPa{x)^ a {x))dx = 6.00. We consider a small 
number of particles so that the breakdown of BCS the- 
ory can be demonstrated clearly in Sec. |III| The small 
particle number means that there are oscillations in the 
Hartree potential in Fig. [TJa) [7j. In the ground state 
the particles occupy the lowest available energy levels as 
shown in Fig. [TJd) [see Eqs. (|7| and 

We consider excitations \v) that are single-particle- like 
in the asymptotic limit but are phase shifted from plane- 
wave scattering states due to the presence of the external 
trapping potential and the trapped supcrfluid gas. For- 
mally, the operators Q\ are constructed from continuum 
one-particle excitations and bound multi-particle excita- 
tions. The multi-particle excitations in general consist of 
all possible configurations of the particles in the Hartree- 
Fock energy levels. We include only the leading-order 
terms, i.e., 

q r,s,t 

where E„ > 0, and E r ,E s ,E t < [8]. The first term 
in Eq. d9| describes continuum one-particle excitations 
[e.g., see Fig. [TJb)]. The second term introduces bound 
two-particle one-hole excitations where the spin of the 
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hole differs from the spin of the projectile particle [e.g., 
the bound excitation with (r, s, t) = (7, 8, 3) is illustrated 
in Fig. life)]. We require E r ,E s > E F and E t < E F as 
these are the only non-zero contributions when acting on 
the Hartree-Fock ground state 

The coefficients C q and B" st in Eq. (|9j) are determined 
so that Eq. is valid, with |G> -> |HF). Taking X^ = 



and — al^a/ s iOti yields 



EgCg 



qtsr 



(10) 



and 



(E r + E s -E t )B» st + J2C"V r 



rstq 



+ ^2 Br's'tVrss'r' ~ ^ st ,V rt > tr > — E v B^ st , 

r' ,s' r',t' 

respectively, where 



(ii) 



v; 



nrnm'n' 



9 I K( x )'Pm( x ) ( t , m'(x)(t>n'(x)dx. (12) 
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FIG. 2: (a) Scattering phase shifts and (b) transmission 
probability for a trapped Fermi gas, calculated using the 
configuration-interaction method. The curves in (a) corre- 
spond to (thick) the even phase shift S c (k) and (thin) the odd 
phase shift S (k). The vertical lines in (b) ind icate the un- 
coupled resonance wavevectors k r3 [see Eq. ( 16 1] . Parameters 
are g = -9.55wE w , V = 127.32£ m , and E F = -19.10E W . 



Equations (10) and (111 define a Hermitian eigenvalue 
problem for the coefficients C q and B^ st . The one- 
particle continuum excitations (C^) are coupled to the 
bound two-particle one-hole excitations (B" st ) by the 
full Hamiltonian ([I]) giving rise to the off-diagonal in- 
teraction terms V q t sr - By construction, the three-particle 
wavefunctions associated with the coefficients B" st vanish 
asymptotically and, in the limit x — > ±oo, the excitations 
\v) are described by the single-particle wavefunctions 



(13) 



To investigate the scattering properties quantitatively, 
the Hartree-Fock ground state is computed on a grid of 
spatial extent 2L. We invoke the Neumann boundary 
conditions [d<p n (x)/dx] x= ±L — so that the Hartree- 
Fock wavefunctions provide an orthonormal basis. Solv- 
ing the eigenvalue problem defined by Eqs. ( 10 1 and (111, 



we determine the single-particle wavefunctions $> v (x) for 
particular discrete values of the energy eigenvalue E v . 
The scattering information at any energy of interest is 
then extracted using the R-matrix method [9], i.e., we 
match Q v {x) to the analytic asymptotic forms 



lim 4>£(x) oc cos(fcx =p SJk)) 

x — >±oo 

lim $°(x) cx sin{kx =p S (k)), 

x — >±oo 



(14) 



where E v — h 2 k 2 /2M, and we use the Hartree-Fock basis 
to reconstruct the asymptotic solutions, and the scatter- 
ing phase shifts S c (k) and S (k), for any k. 

The even and the odd phase shifts for the trapped 
Fermi gas considered in Fig. [I] are shown in Fig.|2^a). We 
observe a variation in the background phase shifts S^%(k) 



due to the effective external potential V oxt (x)+W(x) [see 
Eq. QGJ] . We also observe twenty-four resonance features 
that are characterized by a jump of ~ 7r in cither the 
even or the odd phase shift. 

Transmission and reflection coefficients for the system 
can be calculated from the phase shifts [TO] HI]. For a 
projectile particle incident from x = — oo (or equivalcntly 
from x = oo), the transmission probability (normalized 
to one) is given by |T(fc)| 2 = cos 2 (<5 e (fc) - S (k)). The 
transmission probability is shown in Fig. [2jb). As ex- 
pected, \T(k)\ 2 -> in the limit k -> 0, and |T(fc)| 2 -> 1 
as A; — > oo. The smoothness of the Gaussian exter- 
nal potential means that quantum reflections at the 
trap edges do not play a significant role. Consequently, 
<5j? s (fc) rs <5o S (fc) and the background transmission profile 
approaches unity without the characteristic oscillations 
observed in the case of a square well potential. In addi- 
tion to the background transmission profile, we observe 
twenty- four scattering resonances for which the transmis- 
sion falls to zero. 

The observed scattering resonances are due to interfer- 
ence between the bound two-particle one-hole excitations 
and the one-particle excitation continuum. The coupling 
of the two branches of the excitation spectrum is pro- 
vided by the inter-atomic pairing, i.e., the matrix ele- 



ments V q tsr in Eqs. (10 1 and (111. For each resonance, 



the coupling V q t sr is only non-zero for either the even 
or the odd continuum wavefunctions, depending on the 
parity of the (r,s,t) bound wavefunction [see Eq. ( 12 1] . 



Therefore, a particular (r,s,t) configuration modifies ei- 
ther the even or the odd continuum states, but not both. 

The effect on the continuum of an embedded discrete 
state has been described by Fano [IS]. Following that 
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treatment, the resonant phase shift S^ (k) = S c ^ D (k) 
#eo(fc) near each resonance is approximated by 



tan<5 c R (/c) 



r(*)/2 



where the resonance width is T(k) = 2ML\Vf~t s 
and the resonance energy is determined by 



F t = (hk*„) 2 
2M 



— E r + E s — Et + V rssr — V r i 



(15) 



2 /h 2 k 



(16) 



and x(k) = J2 q \V qtsr \ 2 / (E„ - E q ). Taking V qtsr to be 
approximately independent of q, we find that x(k) — 0. 
The resonance energies are well approximated by E v — 
E* s , as shown by the vertical lines in Fig. |2jb). 

The energy E t rs [see Eq. ( 16 l] can be interpreted as the 



energy required to create the bound two-particle one-hole 
excitation (r, s, t), i.e., to create a spin-up particle in level 
r and excite a spin-down particle from level t to s. The 
last two terms in Eq. (16) account for the accompanying 
change in the pair-wise interaction energy. When the 
energy of the incoming projectile particle matches £y s , 
the bound excitation (r, s,t) is an allowed intermediate 
state for the scattering process [see Figs. [l](b) and [TJc) 
for (r,s,t) — (7,8,3)]. The intermediate state is made 
accessible by the two-body coupling term V q t sr in Eqs. 



(10 1 and (11). 




FIG. 3: (a) The BCS ground state of a trapped Fermi gas. 
The curves correspond to (thick) the Hartree potential U(x), 
(dotted) the pair potential A(x), (thin) the Gaussian exter- 
nal potential V c ^t{x), and (dashed) the combined potential 
V(x) = Kxt(a;) + U(x). (b) A scattering particle excitation, 
and (c) a bound hole excitation, where the energy level oc- 
cupation of the (•) spin up and (o) spin down particles is 
indicated schematically. The arrows indicate the particle re- 
arrangement between the configurations in panels (b) and (c) . 
(d) The ground state occupation of the bo und states in the 
combined potential V(x) [see Eqs. (201 and ( 21 (] is indicated 
by the horizontal component of the markers x . In all panels, 
the horizontal lines indicate (dashed) the chemical potential 
jj, and (solid) the bound state energies E n = Eh < in the 
combined potential V(x). Parameters are g = —9.55wE w , 
V = 127.32B™, and fi = -19.10E W . 



III. MEAN-FIELD APPROACH 

The BCS mean-field treatment of pairing and super- 
fluidity in fcrmionic systems has been extremely success- 
ful. It accurately describes a wide range of systems in 
which many-body pair correlations arc important. The 
tractability of the method is a significant advantage. 
However, it is often difficult to conclusively evaluate the 
validity of the theory. The scattering problem considered 
in Sec. |n]is particularly useful for a discussion of the va- 
lidity of BCS theory because (i) the scattering resonances 
are extremely sensitive to the ground state properties of 
the trapped gas, and (ii) a physical interpretation of the 
resonances is possible allowing the mean-field approach 
to be investigated in detail. 

In the BCS mean-field treatment, we approximate 
Hamiltonian ([TJ by the Hartree-Fock-Bogoliubov Hamil- 
tonian 



H, 



HFB 



/ $,(a:) [Hsp(x) + U{x) - n\ i) a {x)dx 

a J 



H.c. 



(17) 



where the Hartree potential is U(x) = g(tpl l (x)ip a (x)) 
and the pair potential is A(x) = — g(tpi (x)ipi(x)} |13j . 
The two spin states are equally populated so U(x) is spin 
independent. We choose the chemical potential /i < 



so that all of the particles are bound. The Bogoliubov 
transformation , 



n 



(18) 



diagonalizes Hamiltonian (17 1, where the quasiparticle 



amplitudes u n (x) and v n (x) solve the Bogoliubov-de 
Gennes equations [H] 



(19) 



C{x) 


A(x) - 




u n {x) 




u n (x) 


A* (*) 


-L{x) _ 




_ v n {x) _ 


= e« 


v n (x) 



In Eq. (19 1, L(x) = Hsp(x) + U(x) — /i and e„ are the 
quasiparticle energies (taking e„ > 0). The quasiparti- 
cle operators 7^ and 7„ Q obey fermionic commutation 
relations and the quasiparticle modes are populated ac- 
cording to the Fermi distribution function, i.e., at zero 
temperature, (% a %^) = 5 nn >5 a p [Hj. 

Figure [3][ a) shows the self-consistent BCS mean fields 
for a trapped Fermi gas. We have used the same parame- 
ters as in Figs. [T] and [2] However, in the mean-field treat- 
ment, the ground state is not a particle number eigen- 
state [2] . For the parameters used in Fig. [3] the average 
number of particles in each spin state is (N a ) — 5.97 and 
the number variance is (N^) ~ (-^q) 2 = 0.20. 
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To understand the effect of the off-diagonal coupling 
in the Bogoliubov-de Gennes equations, we first consider 
Eq. (19 1 with A(x) — 0. We retain the Hartree potential 
U(x) from the finite A(x) self-consistent solution. In this 
case, Eq. ( 19 1 reduces to the Schrodingcr equation 



~2Mdx 2 



V{x) 



1p n {x) = E n 1p n (x), 



(20) 



where V(x) — V ex t(x) + U(x) is the effective confining 
potential for the many-body system [T5] . Equation ( 20 ) 
has discrete bound states ipb{x) with energy Eb < 0, and 
an excitation continuum of both even and odd scattering 
states ip e k ,0 {x) with energy E k = h 2 k 2 /2M > 0. In the 
particle-hole picture, the excitation spectrum has a par- 
ticle branch [u°(x) = tp n (x)] with energy e° = E n — pt 
(for E n > n) and a hole branch [u°(x) = — t/j n (x)] with 
energy e„ = —E n + /j, (for E n < /i). The superscript zero 
indicates that we are solving Eq. (19 1 with A(x) = 0. 

Taking A(x) to be finite introduces coupling between 
the particle and hole branches of the excitation spec- 
trum. The quasiparticle excitations then become simul- 
taneously particle- and hole-like to reflect the fact that 
the pairing interactions can excite particles to energy lev- 
els lying above the chemical potential. Figure [3jd) shows 
the BCS ground state average particle occupations of the 
uncoupled bound states ipb(%), i.e., 



n h = 



E 



ipb(x)v n (x)dx 



(21) 



The lowest bound levels are fully occupied but there 
is a redistribution of particles near the chemical poten- 
tial, compared to the Hartree-Fock ground state [see Fig. 

Bd)]. 

To investigate the scattering properties of the sys- 
tem, we compute the even and the odd solutions of the 
Bogoliubov-de Gennes equations (19 1 subject to Neu- 



mann boundary conditions. We then match the asymp- 
totic behavior of the quasiparticle amplitudes to their 
analytic forms. The particle-like amplitudes lying in the 
continuum have the form 



lim u c k {x) oc cos(fcx =F <5 c (fc)) 

x — >±oo 

lim u?(x) oc sin(kx =p S Q (k)), 

x — *±oo 



(22) 



where e k = h 2 k 2 /2M — fi. In the asymptotic limit, the 
corresponding hole-like amplitudes v k °(x) tend exponen- 
tially to zero. Again we use the R-matrix method [5] to 
determine the scattering phase shifts <5 c (fc) and S (k) for 
any k. 

The even and the odd phase shifts for the parameters 
used in Fig. [3] are shown in Fig.|4^a). The phase shifts are 
independent of the spin of the projectile particle because 
the two spin states are equally populated. We observe 
a background variation in the p hase shifts due to the 
effective potential V(x) [see Eq. ( 20 1] . We also observe 
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FIG. 4: (a) Scattering phase shifts and (b) transmission prob- 
ability for a trapped Fermi gas, calculated using the mean- 
field method. The curves in (a) correspond to (thick) the 
even phase shift 5 e (k) and (thin) the odd phase shift <5 (fc). 
The vertical lines in (b) indicate the uncoupled resonance 
wavevectors [see Eq. (23 (]. Parameters are g — —9.55wE w , 
V = 127.32£„, and fi = -19.10E W . 



and the odd phase shifts and become narrower higher in 
the continuum. The transmission probability is shown in 
Fig. |3b). 

The observed scattering resonances are possible be- 
cause, in the uncoupled system, a scattering particle ex- 
citation [e.g., see Fig. ^b)] can have the same quasi- 
particle energy as a bound hole excitation [e.g., see Fig. 
|3jc) for b = 3]. In particular, a particle scattering state 
with quasiparticle energy e° = E k — [i is degenerate 
with a bound hole excitation with quasiparticle energy 
e£ = —E^, + /i when e° = e°. In the ordinary particle 
picture this can be rewritten as E k = E£, where 



b\2 



E 



(fifc r b ) 
2M 



-E b + 2\i. 



(23) 



If E° > 0, the uncoupled bound state is embedded in the 
particle excitation continuum and, when A(x) is finite, 
this gives rise to a Fano-type scattering resonance. Note 
that the b = 6 bound state in Fig.[3]does not give rise to 
a resonance because it does not lie sufficiently low in the 
trap, i.e., Eq > 2/i. 

Using standard techniques |16j . we find that near a 
resonance the resonant phase shift is well described by 
the Fano profile 



five resonance features that occur alternately in the even 



tan5 e R (fc) = 
The resonance width is 

r eA k) 2M - 



r e , (fc)/2 



i-4-xeAky 



(24) 



h 2 w k 



O (k)cos(6 h e z(k)-S b o z(k)), (25) 
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where Q ey0 (k) = J ipb(x)A(x)^pl'°(x)dx [T7j and the 
Wronskian W k = ip k (x)dip k (x)/dx — ip%(x)dip^,(x) j dx is 
constant because there is no first order derivative in 



Eq. (20 1 [IS]. For each resonance, the coupling matrix 
clement Q c ,o(k) is only non-zero for either the even or 
the odd continuum wavefunctions, depending on the par- 
ity of the bound hole excitation. The smoothness of the 
Gaussian external potential means that <5 bg (fc) w <5 bs (/c) 
and, therefore, cos(<5 bg (fc) - <5 bg (fc)) « 1. 

The resonance energy is determined by solving e° = 
e b + Xc.o(fc), where 

Xc , ( k ) = 6(fc)±^-Q c 2 (fc)sin( ( 5 e b s(fc)- ( 5 b s(fc)), (26) 

and 9(fc) = f ip h (x)A(x)G k (x, s)A(s)i(j h (s)dxds. The 
upper (lower) sign in Eq. ( 26 1 applies if the uncoupled 
bound hole excitation is even (odd), but sin((5 bg (&;) — 
<5 bg (fc)) « and the first term in Eq. (26 1 dominates. 



The Green's function of Eq. (20 1 is 



G k (x,s) 



M 

Ww k 



(27) 

-+ 0, x(fc) = and the resonances 
the quasiparticle energy of the 



In the limit A(x) 
occur for ejj- = e°, i.e. 
bound hole excitation matches the quasiparticle energy 
of the scattering particle excitation. The resonance en- 
ergy is well approximated by E k = E^ [see Eq. (23)], as 
indicated by the vertical lines in Fig.[3Jb). 

The energy can be interpreted as the energy re- 
quired to excite a bound particle to the chemical poten- 
tial and create a second particle at the chemical poten- 
tial with opposite spin. A scattering resonance occurs 
near this energy if there is coupling between the scatter- 
ing state and the intermediate state, where the projectile 
particle and the bound particle of opposite spin form a 
pair at the chemical potential. In general, this interme- 
diate state is forbidden because the chemical potential 
is not an energy eigenvalue of Eq. (20 1. However, in 



the mean- field theory, the pair potential A(x) facilitates 
pair creation and destruction as if there is a source/sink 
of atom pairs at the chemical potential, i.e., there is an 
effective pair condensate at /i. This allows the projec- 
tile particle and a bound particle of opposite spin to be 
simultaneously removed from the system and the inter- 
mediate state for the scattering process is a bound hole 
excitation with two fewer particles than the scattering 
state [see Figs. [3jb) and|3jc)]. In a number conserving 
treatment, coupling to this intermediate state would be 
forbidden because the Hamiltonian does not couple 
states with different numbers of particles. 

We conclude that, in this case, the scattering reso- 
nances predicted by the mean-field approach are spuri- 
ous. In particular, the mean-field theory does not con- 
serve particle number and this allows for coupling be- 
tween sectors of the Hilbert space corresponding to dif- 
ferent numbers of particles. We discuss the validity of the 
BCS mean-field approach in more detail in the following 
section. 



IV. VALIDITY OF MEAN-FIELD THEORY 

The conventional application of BCS mean- field theory 
provides a steady state ansatz for the quantum system. 
This can be used to compute the expectation values of 
a variety of physical observables that conserve the to- 
tal number of particles. In their seminal paper, Bardeen, 
Cooper, and Schrieffer state explicitly that they only '. . . 
for the moment relax the requirement that the wave func- 
tion describes a system with a fixed number of parti- 
cles . . . ' [19]. They go on to explain how the expecta- 
tion value, of the number conserving Hamiltonian, will 
yield approximately the same result when evaluated ei- 
ther with the BCS state or with the projection of that 
state on any eigenstate of the total number of particles, 
so long as it has a particle number eigenvalue close to the 
(large) mean particle number in the BCS state. 

Our use of BCS theory is less well justified because, 
rather than investigating the full number conserving 
Hamiltonian ([lj within the BCS ansatz, we address the 
dynamics governed by the approximate number non- 
conserving Hamiltonian (17 1. Higher order methods 



based on the number conserving Hamiltonian have been 
applied to scattering from nuclei (e.g., |20j). However, in 
the case of our one-dimensional scattering problem, the 
spurious scattering resonances described in Sec. HI arc 
retained even in such a treatment. 

It is well known that the BCS mean-field approach is 
particularly prone to giving inaccurate results when ap- 
plied to small systems. So, it is perhaps not surprising 
that the BCS treatment predicted spurious resonances in 
Sec. |III| However, by illustrating in detail the nature of 
this breakdown of the theory, we can now present argu- 
ments for why the BCS method can be expected to give 
accurate predictions in higher dimensions and for larger 
systems. 

One of the difficulties with using the BCS mean- 
field approach to describe the one-dimensional scattering 
problem is that the chemical potential is not, in general, 
an allowed energy level in the system. However, this 
problem only arises in systems with discrete energy lev- 
els. For example, in two or three dimensions where the 
energy levels can be highly degenerate, or in a homoge- 
neous system where the energy spectrum is continuous, 
the chemical potential is far more likely to lie at, or very 
near to, an available energy level. Furthermore, the high 
degeneracy provides a compelling argument for the valid- 
ity of the mean-field description in the macroscopic limit, 
as discussed below. 

Figure [5](a) represents the excitation spectrum for a 
generic quantum system where the energy levels are 
highly degenerate near the chemical potential. Apply- 
ing the BCS mean-field theory to such a system would 
give an energy level average population rib that varies 
slowly across the quasi-continuum, as indicated in Fig. 
|5jb). Calculating the scattering properties of the sys- 
tem, in the mean-field approach, we would expect to find 
scattering resonances at energies E k ~ £" b for each sufE- 
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FIG. 5: (a) The energy level diagram for a generic quantum 
system for which the energy levels near the chemical potential 
are highly degenerate, (b) The average energy level popula- 
tion rib that could be expected from a BCS mean-field cal- 
culation. In both panels, the dashed horizontal line indicates 
the chemical potential /x. 



ciently low lying energy level (i.e., for E\, < 2/i). To de- 
termine whether this is reasonable, we could alternatively 
consider this system in the spirit of the configuration- 
interaction approach of Sec. [TTJ In this treatment we 
would expect there to be a large number of resonances 
due to the many possible (r, s, t) configurations. In par- 
ticular, the resonances for a particular value of t would 
have similar resonance energies and, taking the contin- 
uum limit for the r and s levels, the many resonances 
would overlap and we would predict essentially the same 
result as in the BCS mean-field treatment, i.e., a single 
resonance for every low lying bound state t. To verify 
this quantitatively, it would be necessary to show that 
the operator a^a^ in Eq. (jij) could be replaced by a 
complex number that was approximately independent of 
r and s in the quasi-continuum. Therefore, Eq. Q could 
be re-interpreted in the spirit of the Bogoliubov transfor- 
mation (18 1. 



To rephrase the argument, consider a group of A/" states 
near the chemical potential that are on average half filled. 
In this case the states near the chemical potential are 
mutually coupled by the many-body pairing and the cor- 
responding eigenstates can be assumed to be essentially 
symmetric in terms of the population of each single par- 
ticle state. The probability of transferring the projectile 
particle, and a low-lying bound particle of opposite spin, 
to any of the half-filled levels is then amplified by a com- 
binatorial factor increasing with AT, similar to the collec- 
tive spontaneous emission (Dicke superradiance) of light 
from a symmetrically excited atomic medium [21 . The 
collective enhancement would be maximum at exactly 
half filling and the intermediate states populating r and 
s levels near /i would yield one strong resonance for every 



level t, in agreement with the BCS mean-field prediction. 
Additional resonances due to unpopulated r and s lev- 
els lying well above the chemical potential [as in Fig. [2] 
do not benefit from the collective enhacement associated 
with the half filling of levels and would be comparatively 
unimportant in the macrosopic limit. This justifies the 
use of the self-consistent mean-field treatment as a sym- 
metry breaking approach in the same spirit that the semi- 
classical approximation with a mean collective dipole is 
applied to the description of optical superradiance. 



V. CONCLUSION 

We have considered single-particle scattering from 
a trapped two-component degenerate Fermi gas when 
the projectile particle is identical to one of the con- 
fined species. Our theoretical treatment is based on a 
configuration-interaction approach and we predict Fano- 
type scattering resonances that are possible because of 
inter-atomic pairing. The scattering resonances are sen- 
sitive to the ground state properties of the trapped Fermi 
gas and we have described the key features of the scat- 
tering resonances quantitatively. 

We have also presented a BCS mean-field approach 
to the scattering problem and have shown that the non- 
number conservation of the theory leads to spurious scat- 
tering resonances. We have described in detail the break- 
down of BCS theory for this case but have argued that in 
macroscopic systems where the energy spectrum is highly 
degenerate, or continuous, the BCS theory may be ex- 
pected to give accurate results. Furthermore, we have 
suggested that the BCS mean-field approximation can 
be interpreted as an effective semi-classical representa- 
tion of superradiance in the system, and this presents an 
interesting avenue for further research. 

The scattering resonances we predict are relevant for a 
range of reflection and transmission experiments in quasi- 
one-dimensional geometries (e.g., HI])- Transport 
studies of different junction interfaces are of particular in- 
terest for developing electronic devices j25j [26j EZ] • For 
example, the energy sensitivity of the scattering reso- 
nances could allow for energy filtering and, in the case of 
spin imbalance, spin filtering may also be possible. Re- 
cent work has also shown that the scattering properties 
of normal-superfluid interfaces have implications when 
considering the thermodynamics of Fermi gases |28j . 
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